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We examine the effects of disorder on striped phases in high-temperature superconductors and 
related materials. In the presence of quenched disorder, pinning by the atomic lattice - which might 
give rise to commensuration effects - is irrelevant for the stripe array on large length scales. As a 
consequence, the stripes have divergent displacement fluctuations and topological defects are present 
at all temperatures. Therefore the positional order of the stripe array is short ranged, with a finite 
correlation length even at zero temperature. Thus lock-in phenomena can exist only as crossovers 
but not as transitions. In addition, this implies the glassy nature of stripes observed in recent 
experiments. 

PACS numbers: 63.70.+h, 71.45. Lr, 74.72.-h 



I. INTRODUCTION 

Duringptbp last decade evidence emerged from tiieoret- 
ical workfrB and from experiments on cuprateaj and 
closely related nickelatesa lill for the existence of striped 
structures within the MO2 (M = Cu, Ni) planes. These 
stripes are highly correlated states of holes which are 
introduced into the_nlanes by doping and which order 
into a unidirectionafij charge-density wave (CDW, wave 
length a) that may be accompanied by a simultaneous 
spin-density wave of period 2a in the sublatticc_mag=. 
netization of the antiferromagnetic metallic spins .E3£3'L3 
Qualitatively, one may think of stripes as parallel strings 
of holes that constitute an antiphase boundary for spin 
order. 

Particular interest in these stripes arises from the pos- 
sible interplayO between these stripes and supercon- 
ductivity. Thereby it is important to distinguish be- 
tween "dynamic" and "static" stripes. While there is 
evidenca-,fchAt superconductivity can coexist with both 
dynamiccrQ'ta and statio^'Eitdpes, static stripes tend to 
suppress superconductivity,E3iEl in contrast to dynamic 
stripes. Therefore the study of the structure and dynam- 
ics of stripes is of principal importance. 

Various phenomenological pictures have been devel- 
oped for the theoretical description of stripes. While 
charge and spin-order are naturally described within a 
Landau theoryjif] the aspect that stripes act as mag- 
netic domaip-.aialls suggests to describe them as string- 
like objects. E3Ej In the ideally ordered case these strings 
form a periodic array. Dynamic fluctuations are gener- 
ated by thermal and quantum effects, whereas potentials 
tend to suppress dynamic fluctuations while they may 
reduce or increase static conformations of the stripes. 

In such phenomenological models, the crystalline 
structure of the underlying atomic lattice (period b) has 
to be taken into account by a periodic potential, which 
tends to increase the positional order of the stripe, ar- 
ray since is can be the source of lock-in effects. Ido On 
the other hand, the spatially inhomogeneous distribution 



of dopants provides a disorder potential for the stripes 
because of the (screened) Coulomb interaction between 
dopants and holes. For low enough temperatures, these 
dopants can be considered as quenched. 

Pinning by the periodic potential is of particular inter- 
est since lock-in effects might explain the special role of 
certain values for the stripe spacing. In the cuprate sys- 
tem La2- a; Sr a ;Cu04, the mismatch S of magnetic Bragg 
satellite peaks has a nonzero value only beyond a thresh- 
old doping x c « 0.04, where S ss x holds. Eventually a 
saturation value 5 = 4 is observed for x > 0.12.E3 Since 
the mismatch S is related to the lattice and CDW pe=. 
riods via S = b/2a, it allows for a natural explanations! 
of the "a; = ^-problem" :E£l Since at S — | the periods 
of the CDW and the Cu spacings have an integer ratio 
p := adb = 4 this saturation could be a commensuration 
cffect.EZI Similarly, the nickelate system La2- 2; Sr 3 ;Ni04+j / 
shows an o m alous thermodynamic behavior at the values 
p = 2, 3c3E3 which seem to be stable over certain respec- 
tive ranges of the holerGoncentration x + 2yn Even more, 
evidence was reportecfij for plateaus of the mismatch as 
a function temperature at rational values 6 = ^ in 
La 2 Ni0 4 .i25- 

On the other hand, the stripe array can also be pinned 
by disorder .EJo Therefore it is important to take a closer 
look at the competitive pinning by the periodic atomic 
lattice and by disorder in order to understand to what 
extent lock-in effects can persist. A first step in this di- 
rection was made by Hasselmann et alx3 who focus on a 
single stripe. However, since a single stripe and a stripe 
array differ in dimensionality one expects qualitatively 
distinct behavior of the response of the system to disor- 
der. 

The purpose of this paper is to determine the effects 
of periodic and disorder potentials on the structural or- 
der of the two-dimensional stripe array on large scales. 
Quantum fluctuations turn out to be irrelevant in the 
presence of disorder, i.e. the system behaves essentially 
as a classical stripe array. We find that a locked state 
with long-ranged positional order can exist only in the 
absence of disorder. Disorder always leads to unlocking 
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and a state with short-ranged positional order where dis- 
locations proliferate. 

The outline of this article is as follows. In Sec. g 
we establish an elastic model for the quantum stripe ar- 
ray with periodic and pinning potentials. The effects of 
these potentials are discussed in Sec 
cal limit 



In Sec. 



Ill in the classi- 



IV| we demonstrate the irrelevance of 
Our conclusions are drawn in Sec. 



quantum fluctuations 
^ where we discuss the relation of our work to previous 
work, the role of topological defects and implications for 
experiments. 



II. THE MODEL 

We describe the stripe system as an array of interacting 
quantum strings. The strings are assumed to be aligned 
in y direction and to have an average spacing a in x di- 
rection. In the following we ignore topological defects 
in the array, the role of which will be discussed in Sec. 
[y|. Then the stripe array can be considered as an elastic 
system. The displacement field u represents a bosonic 
collective mode of the electron system. Its flu ctuatipns 
are governed by a ( "reduced" ) dynamic actionoO'Ej 
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The action is formulated using the imaginary time (we 
set &b = !)• We identify r = (x,y), /i is a mass den- 
sity and the Hamiltonian H has a contribution from the 
elastic energy 



r -(Vu) 2 

2 V ; 



(2) 



with a stiffness constant 7 which includes the line ten- 
sion of the strings as well as thciu. interaction. Besides a 
contributions of entropic nature|3 the main contribution 
to this interaction will stem from the Coulomb interac- 
tion between the stripes. A further contribution arises 
from a crystal field that aligns the stripes in y direction. 
A priori, the stiffness can be anisotropic with an elastic 



energy density cx ^ x (d x u) 2 + 7 a (<9 a 



To simplify the 



analysis, such an anisotropy can be removed by rescaling 
the y coordinate. Then the effective isotropic stiffness 
constant is related to the original anisotropic constants 
through 7 = y/jxly- Note that j x is dependent on the 
stripe spacing a, i.e., on doping. With increasing distance 
a between the stripes 7 will shrink. 

We will examine the coupling of the stripe array to the 
periodic potential U(x) generated by the atomic struc- 
ture as well as to a random potential V(x,y) due to the 
interaction between the holes and the dopants which may 
be considered quenched at low temperatures. The corre- 
sponding energy contributions read 



Hi 



d 2 rp(r)U(x), H v = / d 2 r p(r)V{x,y) (3) 



in terms of the stripe density 

p(r,r) ~ i{ £ e «3»[-«(^)] - d x u{v, r)} (4) 



where Q m := 2-Km/a are reciprocal lattice vectors of the 
stripe array. The elastic and disorder pinning energies 
of the stripe array are similar to those of vortex lines in 
planar type-II superconductors. A recent review of the 



latter system can be found in Ref. 34 



U is assumed to be periodic, U(x) — U(x + b), with a 
period b < a (the modulation along the stripes is negligi- 
ble for our purposes). For simplicity we take U{x) as an 
even function (this restriction is for the simplicity of our 
analysis but not essential for the results) 



U(x) 



^ U n cos(pQ n x) 



(5) 



n>l 



with p > 1. We assume that the random potential 
V(x, y) is Gaussian distributed with zero average and a 
variance 



V(v)V(0) = 



2tt £ 



(6) 



with a correlation length £ and a weight A. 

Subsequently we will establish the global phase dia- 
gram for the total system with a partition sum 



Z = / V[u] 



(7) 



Since the system without pinning provides an important 
reference point we start with a brief discussion of thermal 
and quantum fluctuations of the displacement. There is 
a characteristic length scale £t ■= y/h 2 j/T 2 [i beyond 
which thermal fluctuations dominate over quantum fluc- 
tuations. A related temperature scale T a := yW^ffaFp, 
is defined by the coincidence It = a. In terms of these 
scales, the displacement fluctuations in systems with a 
large size L> It are obtained as 




27T7 It 2a^/7T7/i 



for T > T a , 
for T< T a . 



(8) 



Thus, while the unpinned stripe array is flat (i.e., (u 2 ) is 
finite for L — » 00) at T = 0, it is logarithmically rough 



(i.e. 



InL for L — > 00) at any finite temperature. 
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III. CLASSICAL LIMIT (h = 0) 

For the analysis of the effects of the potentials it is 
convenient to examine the various limiting cases defined 
by the relative strength of thermal fluctuations, quantum 
fluctuations, periodic pinning and disorder pinning. We 
start from the consideration of the classical limit h — > 
acting on the H appearing explicitly in Eq. ([!]) but not on 
possible implicit dependences of other model parameters. 
In this limit temporal fluctuations become negligible and 
one has to examine the system governed by the Hamil- 
tonian 



TL = TLd + TLu + TLv- 



(9) 



In the absence of the potentials U and V thermal fluc- 
tuations lead to an average displacement that diverges 
logarithmically with the system size L [cf. Eq. (^) for 
T a — 0] which means that the stripe structure has only 
quasi-long-range order in the position of the stripes. 



A. Periodic potential only (V = 0) 

To analyze the relevance of a periodic pinning poten- 
tial we focus on commensurabilities of low order with in- 
teger p. In this case the stripe structure can lock into the 
periodic potential at low temperatures while it unlocks 
at large temperatures. The transition between these two 
states is analogous to the roughening transition of crystal 
surfaces. We follow the standard analysis of the rough- 
ening transition (see Ref. [35] and references therein) in 
order to obtain the transition temperature Tr. Combin- 
ing Eqs. ([|), (||) and (g) we find an average potential 
energy 

(Hu) ~ -L 2 ^ £ U n exp ( -p 2 Q 2 n (u 2 )/2) (10a) 

(10b) 



n — ' 



n>l 



for an infinitesimally weak periodic potential. The low- 
est harmonic n = 1 gives the most relevant contribution 
to this energy. The stripes are locked when the average 
potential energy does not vanish in the limit L — > oo, i.e., 
for temperatures below 



2 7 a 2 
irp 2 



(11) 



The transition temperature increases with increasing 
strength of the potential since 7 is renormalized to larger 
values. The effective parameters on large length scales 



L = ae are ( 
flow equations 



scribed by renormalization group (RG) 



dl 



Ui = (2- 



7tp 



a 2 7 



dl 1 



A—^-U 2 , 



(12a) 
(12b) 



with a temperature dependent coefficient A which is of 
order one near the roughening transition. These equa- 
tions lead to the phase diagram sketched in Fig. [j} 




FIG. 1. Schematic representation of the renormaliza- 
tion-group flow near the roughening transition according to 
Eqs. ([[i]). The line of crosses represents fixed points where 
the periodic potential is irrelevant and the stripes are not 
locked. The dashed line is the phase boundary between the 
locked and unlocked phase. Arrows indicate the RG flow. 

The case with nonintegcr p involves the analysis of 
higher order commensurable states or incommensurable 
states. In such cases the periodic potential is less rele- 
vant than in the low order commensurate cases examined 
above. Since we find later on that the periodic potential 
is irrelevant in the presence of disorder for integer p > 2 
this will be true also for noninteger p. For more de- 
tails on commensurate/incommensurate systems the in- 
terested reader is referred to Rcfs. BqBTl 



B. Pinning potential only (U = 0) 

In order to discuss the relevance of disorder pinning, 
we start from the Hamiltonian Jiy as given in Eq. (j^), 
discard rapidly oscillating terms that are irrelevant on 
scales much larger than a and keep only the most rel- 
evant term m = 1 in the sum over harmonics for the 
density (see Ref. |34| for intermediate steps). After aver- 
aging over disorder we find the effective replica pinning 
Hamiltonian 



a/3 J ^ 



'^cos-K(r)-^(r)] 



(13) 



Disorder couples to the d x u term in the density Eq. (Q) 
as a random field and gives rise to the first term in Eq. 
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( |l3| ) with a value a = (Note that strictly speaking 

this term should contain only d x u a d x u" in the unrenor- 
malized Hamiltonian; the form written in Eq. (|l^) an- 
ticipates that renormalization generates a random field 
coupling to both components of the gradient.) Similar to 
the estimate for the roughening temperature above, one 

,ref 

l v 

respect to 7i c i, 



can estimate the relevance of H.™ p from its average with 



A' 



A 



£ i 2 ^exp(-2^(K-^] 2 )/ a 2 ) 



a, /3=1 



-L 2 N(N - 1) 



A 



where we used (u a (v)vr(v)) = I ^ 2^ — In ^. For temper- 
atures above 



Tsr 



7a' 

TT 



(14) 



the average disorder energy vanishes on large scales. Be- 
low, disorder shows to be relevant and its effects have to 
be calculated by renormalization group techniques. Note 
that 



(15) 



This relation becomes an identity if the renormalization 
of 7 due to the presence of the potentials can be ne- 
glected. Then Tsr >Zr for p > ^2. 

Cardy and OstlunoH were the first to derive RG equa-. 
tions near the transition and Villain and FernandezO 
studied the flow of parameters to their large-scale val- 
ues at zero temperature. A concise summary of these 
two approaches is given in Ref. [}4| We combine the flow 
equations for these two temperature ranges by the inter- 
polation 



da 

~dJ 

dA 

~dJ 



a 2 A 2 
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T 2 7 2 a 4 + A( 7 2 a 4 + c 2 A) ' 
2vrT\ . _ c 2 A 2 



7 a 



0* 



7 2 a 4 + c 2 A 



(16a) 
(16b) 



The numbers c\ and c 2 are of order unity and depend 
only weakly on temperature. 7 is not renormalized due 
to a statistical symmetrypl just like a does not feed back 
to A. This holds for the replica Hamiltonian (|l3|) which 
is a good approximation on large length scales. Smaller 
scales will weakly renormalize the stiffness 7 to larger 
values and generate additional irrelevant terms. From 
the flow equations, Tsr = 7a 2 /tt is identified as the tem- 
perature above which A is renormalized to zero. Never- 
theless, disorder is marginal for T > Tsr since a takes a 
finite fixed-point value. Thus, here one has displacement 
fluctuations 
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2tt \ 7 



(17) 



Below the transition temperature (T < Tsr), A flows to 
a finite fixed-point value and da/dl becomes constant. 
a thus asymptotically has a logarithmic dependence on 
the scale L and gives the dominant contribution to the 
fluctuations 



[u(v)-u(OW 



2 1 2 • 
a %ln -. 

a 



(18) 



Slightly below T SR , x oc (1 - T/T SR ) 2 . This squared- 
logarithmic roughness (|l^ ) defines the superrough (SR) 
phase. In this phase thermal fluctuations can still give 
a logarithmic contribution to the correlator (|l8| ) which - 
however - is sub-dominant. 
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FIG. 2. Schematic representation of the renormalization 
group flow describing the superroughening transition accord- 
ing to Eqs. (0). Arrows indicate the RG flow. The dashed 
line is the phase boundary between the superrough and rough 
phase. Crosses represent fixed points (note, however, that a 
flows to infinity for T < Tsr, = 7a 2 /tt). 



C. With complete potential 



In the presence of both a disorder potential and a pe- 
riodic pinning potential it is not immediately clear which 
one will prevail in what part of the phase diagram. At 
first sight it seems to be possible to have a flat phase, a 
ln-rough phase, or a ln 2 -rough phase. We assume that U 
and V are weak such that a renormalization of the stiff- 
ness 7 as well as of Tr and Tsr is negligible. Then one 
always has Tr < Tsr for p > s/2. Thus, at high temper- 
atures T > Tsr the system will be logarithmically rough 
since both U and V are irrelevant. At intermediate tem- 
peratures in the interval Tr < T < Tsr, V is relevant 
while U is irrelevant with respect to thermal fluctuations. 
This suggests that the system is superrough. Although 
we ultimately find this to be true (see below), the argu- 
ment needs to be refined since it is no longer sufficient 
that U is irrelevant with respect to thermal fluctuations. 
Instead, we one has to argue that U is irrelevant at the 
disorder-dominated, superrough fixed point. Eventually, 
for T < Tr one might expect to find a superrough phase 
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if U is weak compared to V and a flat phase if V is weak 
compared to U. However, the following arguments show 
that a flat phase is not stable in the presence of disorder 
and that superroughness should persist for all T < Tsr. 

In order to argue that a flat phase cannot exist for ar- 
bitrarily weak V and that the stripe array is superrough 
for all T < Tsr, we show (i) that weak periodic potentials 
(U <C V) are irrelevant at the superrough fixed point and 
(ii) that arbitrarily weak disorder (V <C U) will roughen 
the stripe array even for T < Tr, i.e. that disorder is 
relevant at the flat-phase fixed point, which implies that 
U is renormalized to zero and the system is superrough 
on large scales. 

Consideration (i). For U the irrelevance of a peri- 
odic potenti al fol lows from an analysis in analogy to the 
one in Sec. 
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We assume that the disorder-induced 
fluctuations of the displacement field are Gaussian with 
a correlator (u 2 ) ~ ia 2 xln 2 (L/a) [which is implied by 
Eq. dl8|)]. This replaces the thermal correlator in Eq. 
(|lOa|). This now leads to 



(Hi 



n * — ' 



— 7r 2 n 2 p 2 x ln( L/a) 



(19) 



and the average periodic pinning energy vanishes for large 
system sizes. 

Consideration (ii). The situation V <C U is more 
subtle. We assume that the stripes are locked into a 
flat state. We neglect thermal fluctuations which would 
renormalize U to smaller values. For simplicity of our 
argument we assume that p — 2 although the argument 
can be generalized to any p > 1. In addition, we retain 
only the most relevant, lowest harmonic of the periodic 
pinning energy, 



— / a r cos u ( v )- 

a J a 



(20) 



In the absence of V all flat states u(r) = nb with some 
integer n would be equivalent. Disorder certainly breaks 
this degeneracy and one might expect the stripe array 
to find a ground state where it(r) fluctuates only weakly 
around na with some particular n, say n — 0. However, 
one can show that the ground state is not given by small 
fluctuations within this particular "valley" n = but 
that solitons (i.e., local areas where u(r) « Snb with a 
shift Sn 0, cf. Fig. ||) are preferred energetically. The 
proliferation of a large number of such solitons implies 
the irrelevance of U and hence the superroughness of the 
stripe array on large scales. 

We now examine a disk-like soliton of Radius R and 
estimate its elastic energy cost and typical gain of pin- 
ning energy in order to decide whether such solitons are 
favorable. From an energetical point of view the cre- 
ation of Sn = 1 solitons is equivalent to the creation of 
a magnetic domain in a random-field Ising model. For 
a strong periodic pinning potential, U\ ^> 7a, the soli- 
ton has a narrow border of width ijj rs sj 07 /U\b and 



of an energy per unit length e w ^U\a^/p. Thus the 
elastic energy cost is proportional to the border length, 
E c \ oc eR. In the area of the soliton the strings are ex- 
posed to a different disorder potential V (we assume the 
disorder correlation length to be small, £ < a) . Then the 
typical energy gain is proportional to the square root of 
the area Ey oc — y/ Ai? 2 /a£. This gain is larger than the 
elastic energy cost for solitons of all sizes for a disorder 
strength beyond a threshold value 



A c oc a£e . 



(21) 



However, for large solitons one has to take into account 
that the soliton border will be roughened by disorder. 
The equivalent roughening of domain walls in the ran- 
dom field Ising model was studied in Refs. As a 
consequence, the border line tension e is renormalized to 
zero on a finite length scale 



a exp 



<e 2 



(22) 



where c is a constant of order unity. The creation of 
solitons of a size R larger than l t is thus energetically 
favorable. Overlapping solitons of unbounded size imply 
a roughness of the stripe array, provided the sum of the 
shifts ^\ Srii at a given position increases with the num- 
ber of solitons (enumerated by the index i) that include 
this site. In principle, the interaction between solitons 
could lead to a compensation of the shifts between pairs 
of solitons. However, the interaction between the soli- 
tons is short-ranged on the scale l\j and it cannot com- 
pete with a disorder energy that discriminates between 
a shift brii — 1 and a Srii = — 1 for each soliton. Al- 
though the pinning energy of these two states is identical 
in the bulk area, the disorder energy of these two solitons 
is different in the border region and leads to an energy 
contribution proportional to yR. Therefore the multiple 
creation of large and overlapping solitons leads to uncor- 
rected contributions to shifts Sn and therefore implies 
the roughness of the stripe array on scales beyond £ e . In 
this sense U is irrelevant on large scales in the presence of 
arbitrarily weak V, although on small scales the stripes 
will be confined to valleys of U. Since U is irrelevant, the 
stripe array will be superrough as in the absence of U. 

Although the case p — 1 appears not to be of physical 
relevance for striped systems, we add as a side remark 
that for p = 1 a flat phase can exist. In this case the 
elastic energy cost oc R of a soliton cannot be compen- 
sated by a disorder energy which no longer has a bulk 
contribution oc R but only a border contribution oc y/R. 

Strictly speaking, we have shown the irrelevance of U 
only for integer values of p. Since noninteger rational 
values of p would correspond to commensurabilities of 
higher order, they are even more susceptible to the de- 
struction of long-range order by U. 

It is interesting to note that the absence of a flat phase 
for p > 1 is peculiar to two dimensions. In dimensions 
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2 < d < 4 a flat ph-ase is stable for disorder weaker than 
a threshold value. ej'E3 Therefore, a stack of planar stripe 
arrays with a finite coupling between the planes will ex- 
hibit also a flat phase for a disorder strength below a 
certain threshold value. 




2R 

FIG. 3. Illustration of a soliton of radius R for p = 2. Full 
lines represent the strings of holes of average spacing a, dashed 
lines are minima of U with spacing b = a/ p. The soliton width 
is l v . 



IV. ADDING QUANTUM FLUCTUATIONS 

Now we finally study how the potentials U and V affect 
the elastic string model taking into account its quantum 
mechanical nature. In order to analyze whether these 
potentials are relevant at all, the scaling arguments used 
in section 



[II can be applied analogously to the dynamic 



action. Since the displacement correlations of the un- 
pinned array are qualitatively different for T = (where 
the array is flat) and T > (where the array is rough) 
these cases will be discussed separately. 



A. At zero temperature 

In this section we focus on the case with disorder but 
without a periodic potential. Before we turn to the anal- 
ysis of this case, it is worthwhile to point out that this 
system is essentially a two-dimensional "Bose glass." Be- 
cause of an analogy between a <i-dimensional bosonic 
problem at zero temperature and a classical (d + 1)- 
dimensional problem a finite temperature, c3 supercon- 
ducting vortices in the presence of columnar pinning cen- 
ters provide another equivalent system that has received 
a lot of recent attention.E^O The problem at hand, a 
bosonic one-component displacement field at T — 0, is 
equivalent to a stack of classical elastic layers which was 
studied by Balentso According to his analysis - which is 
restricted to dimensions close to d = 4 - the displacement 
field is logarithmically rough. We now focus on two di- 



mensions because characteristic modifications leading to 
superroughness are to be expected. 

The analysis of the disorder pinning is most convenient 
using the replicated action 

S™ P ~Y. J dTdT'd 2 r j-^W*(r,T)Vi/(r,r') 



A 



A[u a (r,T)-u {r, t') 



where 



A(u) 



2tt 
cos — u 
a 



(23) 



(24) 



if we retain only the lowest harmonic of the stripe density 
as most relevant term as in Eq. (|l3|). 

For a scaling analysis of this action contribution, we 
consider a rescaling of space, time, displacement and ac- 
tion quantum according to 



e l r h 



e zl n, 



e (l u h 



(25) 



with a dynamical exponent z, roughness exponent £ and 
action scaling exponent r;. In order to analyze the rel- 
evance of quantum fluctuation later on we allow for a 
rescaling of H. We note that due to a statistical tilt sym- 
metry the flow equation of 7 consists only of the scaling 
parto 



d I 



v + 2Ch 



which implies 



77 = z + 2( 



(26) 



(27) 



at any possible fixed point. If 77 > quantum fluctu- 
ations are irrelevant on large scales and the fixed point 
can be called a "classical" (h^ — 0) fixed point in anal- 
ogy to the. irrelevance of thermal fluctuations for classical 
systems. CJ 

To establish the relevance of pinning we note that the 
action of the unpinned system S e i is invariant under a 
rescaling with z = 1, 77 = and ^ = — |. According to 
Eq. (||) the displacement fluctuations are finite for T = 
(this is reflected by the negative value of C) which implies 
that (cos ^-[u a (r, r) — u@(r, r')]) is finite since in the un- 

h/ {ay/TryJi) for 
oc 



pinned system ([w"(r, r) — u /3 (r, r')] 2 ) 



D 4i 



00 and thus disorder is strongly relevant, (Sy P 



Disorder can be taken into account on a crude level 
(the "random field" approximation) by retaining in ac- 
tion (23) only the harmonic parts bilinear in u. Then 
this action part is found to be scale invariant together 
with 5 e i for z = 1, £ = 1, and 77 = 3. To gain qual- 
itative insight into the spatio-temporal correlations we 
calculate the displacement fluctuations in this random- 
field approximation. We find 



([ U «(r,r)-^(0,0)P) = ^(r,r)+C 2 (r) (28a) 



G 



with a contribution C\ from quantum fluctuations and a 
disorder contribution C2. Cf^(r, r) vanishes for a = (3 
and r = and r = and takes a finite value 



Cf{r, t) 



(28b) 



for (»//Jorr>fflorT> 
roughness 

G T 

C 2 (r)~ -ln- 



'///7a. C*2(r) is rough; its 
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(28c) 



is dramatically overestimated in this approximation. 
From this correlation function we recognize that while 
disorder roughens the displacement in spatial directions 
(consider large r for r = 0), it preserves the flatness in 
temporal directions (consider large r for r — 0). These 
qualitative properties should hold even after renormaliza- 
tion effects due to the anharmonic terms in action 
are taken into account. 

Although a systematic renormalization group analysis 
is very intricate and beyond the scope of this article, we 
present arguments in favor of rj > at the true fixed 
point. First, we use the fact that a fixed point with a 
periodic correlator A(u) can exist only for £ = 0. Thus 



Tj = Z 



(29) 



from Eq. (27) and it is sufficient to show z > for the 
irrelevance of quantum fluctuations. It is natural to as- 
sume that the dynamic exponent is positive, if not diverg- 
ing on large scales as is typical of glassy systems where 
the dynamics is governed by tunneling through divergent 
barriers. _ 

The difficulty^ to calculate z is related to the fact that 
A (it) may become nonanalytic near u — 0. While such 
a nonanalyticity follows from afunctional renormaliza- 
tion group analysis near d — 4,ca it is not clear whether 
such a nonanalyticity is present in d = 2. According 
to the analysis summarized in Sec. 1IIB, such a non- 



analyticity is absent at finite temperature near Tsr. At 
zero temperature, the strong coupling analysis of Villain 
and FernandezEa misses such a possible nonanalyticity. 
Therefore we consider both possibilities. 

We first assume analyticity of the function A(u). In 
this case the "dynamic stiffness" will be renormalized ac- 
cording to 



-/i = (2-* + 2C-ry)M + /yA]. 



(30) 



The functional /3 M [A] represents the vertex corrections 
arising from the Sy P . These corrections are expected 
to be positive since disorder pins the stripes at minima 
of V, thereby confining the temporal fluctuations which 
amounts to an increase of the renormalized \l. If /3 M [A] 
is finite at the fixed point, then 77 = z = 1 + /3 M [A]/(2/i*) 



(//* denotes the fixed-point value of fx). Then 77 > and 
the fixed point will be classical. 

In case the fixed-poinLcorrelator is nonanalytic the 
functional /3 M [A] diverges.E3 This signals a qualitative in- 
crease of the dynamic stiffnesSy-which should be described 
by a kinetic term of the formed 



1 f h/T f 
Skin = ft J dr J d 2 r v\d T u\ 



(31) 



which is more relevant than the original kinetic term in 
Eq. (|i]) . The coefficient v would flow according to 



-v= (2 + C-V> + PA&\ 



(32) 



with another functional /3„[A] > 0. This would imply an 
even larger dynamical exponent z = r\ = 2 + j3 u [A]/v* 
and even stronger irrelevancy of quantum fluctuations. 

Thus, in any case r\ > and the system flows to the 
classical fixed point value for V 7^ = U. Quantum ef- 
fects will result only in a finite renormalization of the 
parameters in the classical system. The most impor- 
tant renormalization effect concerns an increase of the 
dynamic stiffness with a possible generation of v. Al- 
though there is no way of handling a kinetic action of 
the form (|3l|), we expect quantum fluctuations on small 
scales to induce a flat but finite quantum contribution C\ 
to the displacement correlation. The classical contribu- 
tion C2 will be renormalized as in the absence of quan- 
tum fluctuatio ns [t he proper correlation can be obtained 
from equation ( |28q ) by inserting the scale-dependent val- 
ues of A and a as obtained from the flow equations ( |l6| ) 
without rescaling; in the superrough phase a cx In r and 
A cx (In r)/r 2 }. 

Since the presence of the disorder potential implies the 
irrelevance of quantum fluctuations, they cannot be ex- 
pected to modify th e com petition between U and V as 
was analyzed in Sec. Ill C| . Thus, the quantum array has 
superrough spatial correlations at T = 0. 



B. At finite temperature 

An inspection of the correlator (^) of the unpinned 
system suggests that thermal fluctuations dominate over 
quantum fluctuations on large scales. In fact, quantum 
fluctuations are irrelevant at the classical fixed point also 
for T > 0. 

This can be seen from the action as follows. The clas- 
sical fixed points (with both thermal roughness or su- 
perroughness) are described by £ = 0. The finiteness of 
the time integral implies z — and according to Eq. ( pTj ) 
also r\ = 0. Then the effective "dynamical stiffness" flows 
to infinity according to Eq. ( |30| ) or Eq. (j32]). Thereby 
temporal fluctuations of the displacement are suppressed 
on large scales, on which the system is described by the 
static limit 
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(33) 



Thus quantum fluctuations will lead only to a renor- 
malization of the parameters in the classical description. 
Therefore the scaling arguments that were applied in sec- 
tion III to the Hamiltonian hold also for the dynamic ac- 



tion. Thus the system in disorder will be superrough at 
low temperatures - without or with an additional peri- 
odic potential - while it will be thermally rough at high 
temperatures. 



V. SUMMARY AND DISCUSSION 

So far, we have analyzed the stripe array in the elastic 
approximation, i.e. neglecting dislocations. Before we 
discuss the relevance of dislocations, we summarize our 
results. In general, we found quantum fluctuations to 
be irrelevant in the presence of thermal and/or disorder- 
induced fluctuations, i.e., to renormalize the classical 
elastic model only weakly. 

In the absence of disorder (A = 0), the stripe array 
locks into a commensurate periodic potential below a 
roughening temperature Tr = 2jb 2 /ir — 2^a 2 ji^p 2 . This 
transition temperature is proportional to the stiffness 7 
which is implicitly temperature dependent due to an en- 
tropic contribution 7 en tr = nT/a 2 Ei If there were only 
this entropic contribution, the system would always be 
unlocked (T > Tr) at finite temperatures for p > y2- 
However, a temperature independent contribution to 7, 
which arises from the Coulomb interaction, leads to a 
lock-in transition at a finite temperature for p > \/2, 
where the translational order changes from long ranged 
to quasi- long ranged (with logarithmic roughness). 

The large-scale structure of the stripe array is dramat- 
ically influenced by the presence of disorder (A > 0). If 
there were no periodic potential, the stripe array would 
undergo a superroughening transition at Tsr = 70 2 /ir. 
For T > Tsr the array would be unpinned with ln- 
roughness, while it would be pinned and superrough for 
T < ?sr. The same scenario holds also in the presence 
of the periodic potential, from which the stripe array 
always unlocks (on sufficiently large scales even for arbi- 
trarily weak disorder). Therefore, the array is rough at 
all temperatures for A > 0. However, for weak disorder 
and strong periodic potential the crossover length scale 
from flat to superrough correlations will be exponentially 
large, cf. Eq. (|§ 

Hasselmann et al!Ei previously proposed a phase dia- 
gram for a single stripe with flat and disordered phases. 
To our understanding, the disorder was effectively as- 
sumed to have long-ranged correlations, which allows for 
the existence of a flat phase even for a single stripe. In 
contrast to this we consider disorder with short-ranged 
correlations and find that it always roughens the stripe 
array. Because a single stripe represents an elastic system 
of lower dimensionality than the stripe array, our finding 



implies also the roughness of a single stripe in disorder 
with short-ranged correlations. 

Now we come back to discuss the relevance of topolog- 
ical defects in the stripe array, starting with the simplest 
situation for U = V = 0. At low temperatures the stripe 
array can be considered as a "smectic" with quasi-long- 
range translational order and long-range orientational or- 
der. At a temperature 



T m = 



7a 
~8T 



(34) 



it would melt0 due to a proliferation of dislocations into 
a "nematic" liquid with short-range translational order 
and quasi- long-range orientational order, before the pro- 
liferation of disclinations drives a second transition into 
an isotropic liquid with short-ranged orientational order. 
While this scenario is well known for classical systems 
(for a review see e.g. Ref. [32]), its relevance to doped 
Mott insulators was pointed out by Kivelson, Fradkin, 
and EmeryH 

It is instructive to compare the melting temperature 
to the other characteristic temperatures-,related to the 
potentials. As pointed out by Jose et al.^5 distinct melt- 
ing and lock-in transitions can exist only for p > 4 since 
Tr = (4/p) 2 T m . For p < 4 they will merge to a single 
phase transition. The effect of disorder is quite virulent: 
In the elastic approximation, it makes the system su- 
perrough at temperatures below Tsr = 8T m . However, 
superroughness implies that dislocations become energet- 
ically favorable [this follows from the flow equations (|l^) , 
see. e.g. Ref. |5f|. Therefore, their density will be finite 
even at zero temperature. Nevertheless, for weak disor- 
der and at low temperatures the length scale where free 
dislocation appear can be extremely large because of a U 
that tends to lock the stripes up to exponentially large 
scales. Since melting in the absence of disorder occurs 
at a temperature, T m < Tsr, free dislocations will be 
present at all temperatures. 

The presence of dislocations reduces the quasi-long- 
range order (where the satellite peaks have an algebraic 
singularity due to the ln-roughness; ln 2 -roughness corre- 
sponds to an algebraic singularity with an exponent that 
depends weakly on the wave vector) to short-range or- 
der with a correlation length of the order of the distance 
between (free) dislocations. Thus, the saturation of the 
dislocation density at low temperatures implies also the 
saturation of the correlation length at a maximum value. 
This conclusion is consistent with experimental observa- 
tions in the cupratesE3 (here it was observed in magnetic 
ordering, which can be destroyed due to a coupling to 
disorder even without a-distortion of the stripe array) as 
well as in the nickelatesESEl (here the correlation length 
is explicitly that of charge order). 

Note that the true correlation length of charge order £ c 
(which we identify with the distance between free dislo- 
cations) can be related to the measured width of peaks in 
the structure function only if the wave-vector resolution 
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is much smaller than l/£ c . Thus, a system without dislo- 
cations (with In- or ln 2 -roughness) has an infinite correla- 
tion length and the apparent correlation length deduced, 
from experiments would be resolution limited. ZachaiO 
recently proposed an explanation of the observed appar- 
ent correlation lengths in terms of chaotic fluctuations of 
the distance between neighboring stripes (see also Ref. 
|30| ), excluding explicitly a key role of dislocations. How- 
ever, this argument is based on the assumption of non- 
integer p and cannot account for the finite . c a rr . elation 
length observed for integer p in the nickelatesO£!rLJ and 
the cupratesEa It also cannot explain why the transverse 
correlation length (along-,the stripes) - which apparently 
is not resolution limitedEa - can be finite. Even more, the 
fact that longitudinal and teansyerse correlation lengths 
are roughly of the same sizeEaH is consistent with their 
relation to the dislocation density. 

As both the roughening transition is washed out by 
disorder and the superroughening transition is washed 
out by the presence of dislocations, no sharp transi- 
tion will exist in the thermodynamic limit. Neverthe- 
less, crossover phenomena may be observable. For weak 
disorder and/or strong lattice potential the stripe ar- 
ray appears to be locked to the lattice potential up to 
very large length scales. It "unlocks" also on finite scales 
with increasing temperature or decreasing lattice poten- 
tial. An apparent lock-in transition was observedfl in 
L^-z-yNdySr^CuC^ where the strength of the lattice 
potential corresponds to the Nd concentration. Since 
2s r S> T m the density of free dislocations will be high 
near Tsr so that reminiscences on finite scales of the su- 
perroughening transition are unlikely to survive. 

Since the array has short-range order at all temper- 
atures for A > we expect the absence of commensu- 
rate/incommensurate transitions (identified from fluctu- 
ations on asymptotically large length scales). Neverthe- 
less, the experimental observations of anomalies at par- 
ticular values of p can be related to the effects of U on 
finite length scales. 

In view of our conclusion about the absence of com- 
mensuration effects we comment on the S = | problem, 
i.e., the observation that S — l/2p apparently saturates 
at S = | near a doping x = 0.125 of the cuprates. If thia 
were a commensuration effect as suggested previouslyJj'O 
one would expect to observe the value 6 = g around 
x = |, i.e. above (x > S) and below (x < S) the match- 
ing x = S. To the best of our knowledge (cf. the data 
collected in Fig. 7 of Ref. |25|), there is no evidence fox 
data with x < S. As pointed out by Yamada et al.,tB 
the saturation of 5 coincides with a saturation of the ef- 
fective hole concentration in the CuC>2 planes beyond a 
certain doping level x « 0.12. Thus it is conceivable that 
other mechanisms limit the effective hole concentration 
and lead to a plateau in d. 

In conclusion, we have pointed out the relevance of dis- 
order for a stripe array even in states where its period 
is commensurate with the atomic structure. We found 



that on large scales pinning by disorder dominates over 
pinning by the atomic structure. This induces the super- 
roughness of the array and, on sufficiently large scales, 
the presence of free dislocations even at low tempera- 
tures, which explains the saturation of the correlation 
length observed in experiments. 
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